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Detecting Sunyaev-ZeFdovich clusters with Planck: 
I. Construction of all-sky thermal and kinetic SZ-maps 
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ABSTRACT 

All-sky thermal and kinetic Sunyaev-Zel'dovich (SZ) maps are presented for assessing how 
well the Planck mission can find and characterise clusters of galaxies, especially in the pres- 
ence of primary anisotropies of the cosmic microwave background (CMB) and various galac- 
tic and ecliptic foregrounds. The maps have been constructed from numerical simulations of 
structure formation in a standard ACDM cosmology and contain all clusters out to redshifts 
of z — 1.46 with masses exceeding 5 x lO 13 M0//z. By construction, the maps properly ac- 
count for the evolution of cosmic structure, the halo-halo correlation function, the evolving 
mass function, halo substructure and adiabatic gas physics. The velocities in the kinetic map 
correspond to the actual density environment at the cluster positions. We characterise the 
SZ-cluster sample by measuring the distribution of angular sizes, the integrated thermal and 
kinetic Comptonisations, the source counts in the three relevant P/ancfc-channels, and give the 
angular power-spectra of the SZ-sky. While our results are broadly consistent with simple es- 
timates based on scaling relations and spherically symmetric cluster models, some significant 
differences are seen which may affect the number of cluster detectable by Planck. 

Key words: galaxiesxlusters: general, cosmology: cosmic microwave background, methods: 
numerical, space vehicles: Planck 



1 INTRODUCTION 



The Sunvaev-Zel'dovich (SZ) effects ISunvaev & Zerdovichll972L 
ll98(ilBirkinshawlll999l:lReDhaelilll995h have evolved from phvs- 
ical peculiarities to valuable and sound observational tools in cos- 
mology. The thermal SZ-effect arises because photons of the cos- 
mic microwave background (CMB) experience Compton-collisions 
with electrons of the hot plasma inside clusters of galaxies and 
are spectrally redistributed. The amplitude of the modulation of 
the Planckian CMB spectrum is a measure of the cluster electron 
column density and temperature. Alternatively, CMB photons may 
gain energy by elastic Compton collisions with electrons of the 
intra-cluster medium (ICM) due to the peculiar motion of the clus- 
ter relative to the CMB. This so-called kinetic SZ-effect is propor- 
tional to the peculiar velocity weighted electron column density 
and directly measures the cluster's velocity component parallel to 
the line-of-sight relative to the comoving CMB frame. 

The advancement in sensitivity and angular resolution of sub- 
millimeter and microwave receivers have allowed high quality in- 
terferometric imaging of more than fifty clusters of galaxies by 
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ground based telescopes ICarlstrom et alj2002h out to redshifts of 
~ 0.8, despite incomplete coverage of the Fourier plane. Apart from 
its primary scientific objective, namely the cartography of the CMB 
with angular resolutions close to 5', th e upcoming Planck mission 
(Mandolesi et al. 1995; Tauber 2000) will be an unique tool for 
observing clusters of galaxies by their SZ-signature. Planck is ex- 
pected to yield a cluster catalogue that is surpassing the classic op- 
tical Abell catalogues or any existing X-ray catalogue in numbers 
as well as in depth and sky coverage. 

The capability of Planck to detect SZ-clusters has 
been the subject of many recent works, pursuing ana lytical 
jAghanim et alJll997t iDelabrouille et aljEoOa lBartelmanr]|200lt 
Moscardini et alJ 120021) as well a s semi analytical iKav et all 
1200 it iHobson & McLachlarJ l2003t) and numerical approaches 
iHerranz et alj2002HWhitel2003t) . Their consensus is an expected 
total number of a few times 10 4 clusters and the detectability of 
(sufficiently massive) clusters out to redshifts of z ~ 1. The authors 
differ mainly in the expected distribution of the detectable clusters 
in redshift z. Where adressed, the authors remain sceptic about the 
detectability of the kinetic SZ-effect. 

As a result of various approximations made, there are clearly 
limitations in these studies: Concerning the SZ profiles of isolated 
clusters, simplifying assumptions like spherical symmetry, com- 
plete ionisation and isothermality have usually been made. Ana- 
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lytical treatments mostly rely on y8-profiles for modeling the spatial 
variation of the Compton-y parameter. Temperature models mostly 
make use of scaling laws derived from spherical collapse theory 
or are taken from X-ray observations. Naturally, the halo-halo cor- 
relation function is not taken account of, neither do the velocities 
correspond to the actual density environment, they are commonly 
drawn from a (Gaussian) velocity distribution. 

The primary application of the all-sky SZ-maps would lie in 
the assessment of the extent to which cluster substructure and de- 
viations from spherical symmetry, the halo-halo clustering and de- 
viations from the scaling-laws alter the predictions made based on 
analytic methods. 

Additionally, the investigations mentioned above lack the in- 
clusion of galactic foregrounds (f or a com prehensiv e revie w of 
foregrounds concerning Planck, see Bouchet & Gispert 1999), the 
thermal emission from planets and minor celestial bodies of the so- 
lar system, beam patterns and spatially non-uniform instrumental 
noise. In order to quantify the extent to which the galactic and eclip- 
tic foregrounds impede the SZ-observations by Planck, i.e. down to 
which galactic latitudes clusters will be detectable, a detailed sim- 
ulation is necessary. Furthermore, the noise patterns will be highly 
non-uniform due to Planck's scanning strategy. For investigating 
this issue, all-sky maps of the thermal and kinetic SZ-effects are 
essential. 

In our map construction, we take advantage of two numeri- 
cal simulations of cosmic structure formation: The Hubble-volume 
simulation, that provided a well-sampled cluster catalogue covering 
a large volume and secondly, a set of template clusters resulting 
from a gas-dynamical simulation on much smaller scales, allow- 
ing us to extract template clusters. For all clusters of the Hubble- 
volume simulation, a suitable template was chosen and after having 
performed a scaling operation to improve the match it is projected 
onto the celestial sphere at the position requested by the Hubble- 
volume catalogue. By construction, the resulting all-sky SZ-maps 
show halo-halo correlation even on large angular scales, incorpo- 
rate the evolution of the mass function and have the correct size 
distribution. In the kinetic SZ-map, it is ensured that the cluster 
peculiar velocities correspond to the ambient cosmological density 
field. Furthermore, the template clusters do exhibit realistic levels 
of substructure and departures from isothermality, and their ensem- 
ble properties also account for scatter around the idealised scaling 
laws. Therefore, most of the imperfections of traditional approaches 
will be remedied by our map construction process. However, there 
are impediments that could not be disposed of: They include gas 
physics beyond adiabaticity, e.g. radiative cooling and supernova 
feedback, that significantly alter the baryon distribution and tem- 
perature profiles of the ICM and hence the SZ-amplitude, incom- 
plete ionisation, inclusion of filamentary structures and uncollapsed 
objects or diffuse gas. Another process influencing the thermal his- 
tory of the ICM is reionisation, which also had to be excluded. Yet 
another complication are non-thermal particle populations in clus- 
ters of galaxies that give ris e to the relativistic SZ-effect IWrightl 
ll979HEnfflin & Kaisej2o69) . 

The paper is structured as follows: After the definition of 
the basic SZ quantities in Sect. |2| the simulations are outlined in 
Sect.[3] The construction of the maps is described in detail in Sect.|3] 
and the properties of the resulting maps are compiled in Sect.|5] Fi- 
nally, the conclusions are presented in Sect.|6| 



2 SUNYAEV-ZEL'DOVICH DEFINITIONS 

The thermal and kinetic Sunyaev-Zel'dovich effects are the most 
important sources of secondary anisotropies in the CMB. Compton 
interactions of CMB photons with electrons of the ionised ICM 
give rise to these effects and induce surface brightness fluctuations 
in the CMB sky, either because of the thermal motion of the ICM 
electrons (thermal SZ) or because of the bulk motion of the cluster 
itself (kinetic SZ). 

The relative change AT IT in thermodynamic CMB temper- 
ature at position as a function of dimensionless frequency x = 
hv/(k B TcMB) due to the thermal SZ-effect is given by eqn. Q: 

AT I e" + 1 \ 

—(9) = y(fl)U_^-4l with (1) 

y(0) = f dln c (0J)TA6J), (2) 

where the amplitude y of the thermal SZ-effect is commonly known 
as the thermal Comptonisation parameter. It is proportional to the 
line-of-sight integral of the temperature weighted thermal electron 
density (c.f. eqn. 0). m c , c, k B and cr T denote electron mass, speed 
of light, Boltzmann's constant and the Thompson cross section, re- 
spectively. The kinetic SZ-effect arises due to the motion of the 
cluster relative to the CMB rest frame parallel to the line of sight. 
The respective temperature change is given by: 

— (0) = - w (0) with w(0) = — fdl njfi, l)v r (0, 1). (3) 
T c J 

Here, v r is the radial component of the cluster velocity, i.e. the ve- 
locity component parallel to the line-of-sight. The convention is 
such that the CMB temperature is increased, if the cluster is mov- 
ing towards the observer, i.e. if v, < 0. In analogy to y, the quantity 
w is refered to as the kinetic Comptonisation parameter. 



3 SIMULATIONS 

Due to the SZ-clusters being detectable out to very large redshifts, 
due to their clustering properties on very large angular scales, and 
due to the requirement of reducing cosmic variance when simulat- 
ing all-sky observations as will be performed by Planck, there is the 
need for very large simulation boxes, encompassing look-back dis- 
tances to redshifts of order z — 1 which corresponds to comoving 
scales exceeding 2 Gpc. Unfortunately, a simulation incorporating 
dark matter and gas dynamics that covers cosmological scales of 
that size down to cluster scales and possibly resolving cluster sub- 
structure is presently beyond computational feasibility. 

For that reason, a hybrid approach is pursued by combin- 
ing results from tw o simulations: The Hubble-volume simulation 
I Jenkins et al. 2001; Colberg et al. 2000J), and a smaller scal e sim- 
ulation including (adiabatic) gas physics (White et al. 2002). The 
analysis undertaken by Bartelmann 1 2001 ) gives expected mass and 
redshift ranges for detectable thermal SZ-clusters, which are cov- 
ered completely by the all-sky SZ-map presented here. 

The assumed cosmological model is the standard ACDM cos- 
mology, which has recently been supported by observations of the 
WMAP satellite ISpergel et all2003l) . Parameter values have been 
chosen as fi M = 0.3, S1 A = 0.7, H = lOO/jkms" 1 Mpc~' with 
h = 0.7, n B = 0.04, n s = 1 and cr 8 = 0.9. 
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data set 


z range 


n/sr 




shape 


MS 


0.0 sS z i 0.58 


4tt 


564818 


sphere 


NO 


0.58 <z< 1-46 


n/2 


182551 


northern octant shell 


PO 


0.58 <z< 1.46 


n/2 


185209 


southern octant shell 



Table 1. Basic characteristics of the light-cone outputs used for compiling a 
cluster catalogue. The columns denote the label of the data set, the range in 
redshift z, the solid angle Q. covered, the number of objects A^i,, retrieved, 
and the geometrical shape. 



3.1 Hubble-volume simulation 

The Hubble-volume simulation is one of the largest simulations 
of cosmic structure formation carried out to date. The simulation 
domain is a box of comoving side length 3 Gpc//? (for the stan- 
dard ACDM cosmology) and comprises 10' dark matter particles. 
The simulations used were carried out by the Virgo Supercomput- 
ing Consortium using computers based at the Computing Centre of 
the Max-Planck Society in Garching and at the Edinburgh parallel 
Computing Centre. The data are publicly available for download 1 . 

The light-con e output of the Hubble-volume simulation 
jEvrard et al1l2002l) was used for compiling a cluster catalogue. 
This ensures, that the abundance of clusters at any given redshift 
z corresponds to the level of advancement in structure formation up 
to this cosmic epoch. The minimal mass was set to 5 x 10 13 Mq//i, 
which roughly corresponds to the transition mass between a rich 
group of galaxies and a cluster. In order to cover redshifts out to 
the anticipated limit for Planck, light-cone outputs of differing ge- 
ometry were combined: First, a sphere covering the full solid angle 
of 4tt was used for redshift radii of z = to z = 0.58. For red- 
shifts exceeding z = 0.58, the northern and southern octant data 
sets were added. The octant data sets span a solid angle of n/2 and 
were replicated by rotation in order to cover the full sphere. TableQ 
summarises the properties of the different output geometries. In this 
way, a cluster catalogue with cluster mass M, position on the sky 0, 
redshift z and peculiar velocity v T projected onto the line-of-sight 
was compiled, comprising a total number of 2035858 clusters. For 
the sky-map construction, the positions 9 were interpreted as eclip- 
tic coordinates, the default coordinate system for Planck. 

Here, it should be mentioned that the combination of different 
outputs gives rise to boundary discontinuities, at the surface of the 
central sphere as well as on the faces of the octant shells. These dis- 
continuities do not only show up in the spatial halo distribution, but 
also in the velocities of clusters close to simulation box boundaries. 
Furthermore, the cluster catalogues exhibit small completeness de- 
ficiencies close to the edges of the simulation domain. 



3.2 Small scale SPH cluster simulations 

A hydrodynamical simulation of cosmological structure formation 
(White etalJl200l constitutes the basis of the SZ-template map 
construction. The simulation was performed with the GADGET code 2 
I Sorineel et al. 2001) using th e 'entr opy-conserving' formulation 
of SPH l Sorinael & Hernquist 2002). The simulation, first anal- 
ysed in White et al. 1 2002), followed 216 3 dark matter particles as 
well as 216 3 gas particles in a cubical box of comoving side length 
100 Mpc//i with periodic boundary conditions. Purely adiabatic gas 



physics and shock heating were included, but radiative cooling and 
star formation were ignored, which however does not result in sig- 
nificant differences in SZ morphology, as shown bv IWhite et alJ 
(2002), but would impact on the scaling relations as demonstrated 
bv lda Silva et alJ 1200 ll) . We analyse 30 output redshifts ranging 
from z = out to z = 1.458. The comoving spacing along the 
line-of-sight of two subsequent outputs is 100 Mpc//j. Halos were 
identified using a friends-of-friends algorithm with linking length 
b = 0.164, which yields all member particles of cluster-sized 
groups. We then employed a spherical overdensity code to estimate 
the virial mass and radius of each cluster. We computed the mass 
M vir inside a sphere of radius r vir , interior to which the average 
density was 200 times the critical density p cr j t = 3H(z) 2 /(SnG). We 
imposed a lower mass threshold of M vlr > 5 x 10 13 M o //i in or- 
der tomatch the lower mass limit adopted for the Hubble-volume 
cluster catalogue. 



4 SUNYAEV-ZEL'DOVICH MAP CONSTRUCTION 

The construction of the all-sky SZ-map proceeds in three steps: 
First, a set of template cluster maps was derived based on cluster 
data from a gas-dynamical simulation (Sect. l4~TY Then, for each 
of the clusters in the cluster catalogue obtained from the Hubble- 
volume simulation, a suitable hydrodynamical cluster template has 
been selected, scaled in mass and temperature in order to better fit 
the cluster from the Hubble-volume catalogue (Sect. 1431 . and, for 
the kinetic sky map, boosted to the radial peculiar velocity required 
by the Hubble-volume simulation. The last step is the projection 
onto a spherical celestial map (Sect. 1431 . In the subsequent para- 
graph (Sect. l4"4*l . the completeness of the resulting SZ-maps is in- 
vestigated analytically. 



4.1 SZ-template map preparation 

Square maps of the Compton-y parameter of the selected clusters 
were generated by SPH projection of all friends-of-friends iden- 
tified member gas particles onto a Cartesian grid with 128 2 mesh 
points. The (comoving) side length s of the map was adapted to the 
cluster size, such that the comoving resolution g = s/128 of the 
grid is specific to a given map. 

If the particle p at position r p = [x p ,y p ,z^ has a smoothing 
length h p , an SPH electron number density n p , and an SPH electron 
temperature T p , the Compton-y parameter for the pixel at position 
x is given by: 



y(x) 



m e c 2 g 



x+g/2 y+g/2 



x-g/2 y-g/2 



with r ■ 



yj(x p - x) 2 + (y p - yf 



(4) 



(5) 



http: //www.mpa-garching .mpg .de/galform/virgo/hubble 



http : //www . mpa-garching .mpg . de/galform/gadget/index . shtml 



Here, we assume complete ionisation and primordial ele- 
ment composition of the ICM for the determination of electron 
number density and temperature. In this way, we produce pro- 
jections along each of the three coordinate axes. The function 
7C is the spherically symmetr ic cubic spline kernel suggested by 
iMonaehan & Lattanzid ll985h . which was also used in the SPH 
simulation: 

g f 1 - 6u 2 + 6m 3 , < u < 1/2 
<K(u) = - x | 2(1 - uf , 1/2 < u < 1 with u = r/h„. (6) 
U ,u> 1 
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Figure 1. The population of our template clusters in the mass-redshift plane. 
The line separates the sparsely sampled region from the region in which a 
sufficient number of template clusters is available. 



Figure 2. Distribution of the ratios ay (circles, solid line) and ay/ (crosses, 
dashed line) that describes the impact of the mismatch between friends-of- 
friends masses and virial masses on the Comptonisations J/ and 'W. 



The fact that the kernel % has a compact support u e [0 ... 1] 
greatly reduces the computational effort. 

The kinetic maps were treated in complete analogy: Maps of 
the Thomson optical depth r were derived by means of eqn. 



t(x) 



h 3 

<4i 



x+g/2 y+g/2 



fafafa*(ih 



x-g/2 y-g/2 



(7) 



In eqn. 0, the influence of velocity differences inside the 
clusters was omitted. At the stage of projecting the template clus- 
ters onto the spherical map, the r-map obtained is boosted with 
the peculiar line-of-sight velocity u, in units of the speed of light 
c required by the entry in the Hubble- volume catalogue in order to 
yield a Compton-w amplitude. 

Neglecting velocity differences inside the clusters does not se- 
riously affect the measurement of cluster peculiar velocities with 
the kinetic SZ-effect shown by Naeai et al. 1 2003), the scatter in 
the velocity estimates increases only little (50 - 100 km/s) when 
considering a rather narrow beam (A0 = 1 '0 FWHM), while the 
kinetic SZ-amplitude remains an unbiased estimator of the pecu- 
liar velocity. For our application purpose, the situation is even less 
troublesome because of Planck's wide beams ( i 5'.Q (FWHM)). 

In this way, a sample of 1518 individual template clusters was 
obtained, and maps for projections along all three coordinate axes 
were derived, yielding a total of 4554 template maps for each of 
the two SZ-effects. Fig. Q shows the distribution of clusters in the 
mass-redshift plane. Especially at high masses, the smooth growth 
of clusters by accretion can be clearly seen. Sudden jumps to larger 
masses are caused by the merging of low-mass clusters. 

A limitation to our SZ-map construction is immediately ap- 
parent: The 4554 cluster template maps derived from the hydro- 
simulation are not strictly independent, but merely show the same 
clusters at different redshift. Thus, the morphological variety re- 
mains limited, but even though there is of course some variation in 
morphology due to accretion and merging events. This, however, 
may not be a severe restriction, keeping the wide Planck-beams in 
mind, that are unlikely to resolve cluster substructure for a large 
fraction of detectable clusters. In this case, the simulation will pick 
up mismatches in Comptonisation relative to the value expected 



from the spherical collapse model in conjunction with the Press- 
Schechter distribution of halo masses. 



4.2 Cluster selection and scaling relations 

In order to select a template map for projection, the closest template 
cluster in the log(M)-z-plane for a given cluster from the Hubble- 
volume simulation was chosen. For the sparsely sampled region of 
the M-z-plane to the right of the line in Fig.Q a cluster from a pool 
containing the most massive clusters to the right of this line in the 
redshift bin under consideration was drawn. 

The template clusters are scaled in mass, temperature and spa- 
tial extent in order to yield a better match to the cluster from the 
Hubble-volume simulation according to formulae l8H0l The scaling 
is parameterised by the masses of the cluster of the Hubble-volume 
simulation yvf (Hubblc) an( j the template cluster yVf". cmplatc> : 

vir r vir 



M 



9r 



q T 



(template) 



vi: 

n (m. 



, Hubble)^ +f . 3 



(template) 



(M« ubbk) )' 2 + , 



h (m. 



(template) V 2 

I 3 



(8) 



(9) 



(10) 



The parameters r, and t„ i e (1,2,3), describing the scaling 
in size q R and in temperature q T were derived from template data: 
Fits to the virial radius as a function of mass and of the mean tem- 
perature inside the virial sphere as a function of mass were ap- 
plied to the data of simulation outputs binned in five data sets. 
This approach leaves the map construction independent of idealised 
assumptions, like the prediction of cluster temperatures from the 
spherical collapse model, or from electron temperature measure- 
ments deduced from X-ray observations and keeps the weak trend 
of cluster temperature with redshift z as contained in the simula- 
tions. 

Although the scaling has been constructed in order to yield 
the best possible match between the template cluster and the target 
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cluster from the Hubble-volume simulation, there are artifacts in 
irregular systems due to inconsistencies in cluster masses A/£" p 
determined with a friends-of-friends algorithm for identifying clus- 
ter member particles and the virial mass estimates ^ tem ? late) follow- 

r vir 

ing from applying the spherical overdensity code. After the scaling, 
these mismatches may be expressed as: 



and 



en) 

(12) 



FiglH shows the distribution of ratios ay and a<w for the en- 
tire Hubble-volume catalogue. Clearly, one recognises large tails 
towards high values, because clusters have on average to be scaled 
to higher masses. This is due to the fact that the hydro-simulation 
outlined in Sect. 13.21 does not sample the high-mass end of the 
Press-Schechter function satisfatorily, simply because of its small 
volume. Nevertheless, the mean of the distributions is close to one, 
which implies that the mismatches average out for the bulk of clus- 
ters. 



4.3 Projection onto the celestial sphere 

For storing all-sky maps the HE ALPix 3 tesselation of the sphere 
proposed by Gorski et al. 1 2002) has been chosen. In order to sup- 
port structures as small as clusters, the Af sidc parameter has been set 
to 2048, resulting in a total number of 12 = 50331648 pixels. 
The side length of one pixel is then approximately 1(71, which is 
well below the anticipated Planck beam size of 510 in the highest 
frequency channels. 

The scaled cluster maps are projected onto the spherical map 
by means of stereographic projection at the south ecliptic pole of 
the celestial sphere. By dividing the (comoving) position vector 
(x, y) of a given pixel on the template map by the comoving angu- 
lar diameter distance x(z) at redshift z, one obtains the coordinates 
(a,[}) on the tangential plane. Then, the stereographic projection 
formulae yield the (Cartesian) position vector (£, rj, of this point 
projected onto the unit sphere: 



4a 





4 + a 2 + p 1 ' 4 + a 2 + f} 2 ' 4 + a 2 + fi 2 



4/8 



(13) 



In order to assign a Comptonisation amplitude to a given 
HEALPix pixel in the projection process, a solid angle weighted 
average is performed. For close-by clusters, the mesh size of the 
templates converted to angular units is larger than the HEALPix 
pixel scale. For those clusters, the map is refined iteratively by sub- 
division of a pixel into 4 smaller pixels subtending a quarter of the 
original solid angle until the pixel size is well below the HEALPix 
pixel scale. Before projection, the template maps are smoothed 
with a Gaussian kernel with AO = 2'.0, which is comparable to the 
HEALPix pixel scale. In this way, it is avoided that structures are 
destroyed by the combination of multiple template map pixels into 
a single HEALPix pixel. This convolution does not affect the later 
usage for simulations concering Planck: A second successive con- 
volution with the narrowest beam results in an effective smoothing 
of 5.38', which corresponds to a decrease in angular resolution of 
roughly 7.5%. 

Additionally, a rotation of the template map around the e - 
axis about a random angle is performed in order to avoid spurious 
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alignments of clusters. The projected pixels are then transported by 
Euler-rotations of the vector (£, rj, f ) to the position requested by 
the Hubble catalogue. 

4.4 Completeness of the all-sky SZ-maps 

The angular resolution of Planck will not allow to spatially resolve 
low-mass and high-redshift clusters. There will be a Compton-y 
background (yb g >s due to the higher number density of low-mass 
clusters compared to high-mass clusters which overcompensates 
their lower individual SZ-signature. Since ideally any isotropic 
background could be removed, we only have to take into account 
the average background fluctuation level {yl s )g which is described 
by power spectrum statistics. 

This section studies the influence of the background of unre- 
solved SZ-clusters in our all-sky map of SZ-clusters on power spec- 
trum statistics. Our simulation neglects the SZ signal of clusters 
both with masses smaller than 5 x 10 13 Mq/7? and redshifts z > 1.5. 
In principle, these missing clusters could be accounted for by draw- 
ing them from a particular realisation of a Poissonian random field 
such that they obey the correct relative phase correlations, i.e. that 
they exhibit the observed local clustering properties. 

However, there are large uncertainties about the baryon frac- 
tion / B = £2b/£2m of low-mass halos(M hato < 5 x 10 13 Mo/h) es- 
pecially at higher redshifts. Analyses of X-ray observations of 45 
lo cal clusters (z < .18, only 4 of them lie at z > 0.1) carried out 
bv lMohr et all <1999l) and lArnaud & Evrard <1999l) suggest a weak 
trend of the cluster baryon fraction / B with cluster mass M and 
a deviation from the universal value, which may be due to feed- 
back processes like galactic winds that more effectively deplete the 
ICM of baryons in low-mass compared to high-mass clusters. The 
behaviour of fs at high redshifts is very uncertain, among other 
reasons because the relative importance of the different feedback 
processes at high redshift is yet unknown. This uncertainty is also 
reflected in different cooling rates and mechanisms, governing the 
ionisation fraction of the electrons and the resulting SZ flux of a 
particular cluster. In the following, we study the contribution to the 
SZ flux of clusters both with masses smaller than 5 x 10 13 Mq//i 
and redshifts z > 1.5. Although the impact of this cluster popu- 
lation to the (y bl! )-statistics amounts to a significant fraction, this 
population has a negligible contribution to the more relevant (y 2 ,,)- 
statistics which will be shown in the following. The unresolved 
cluster population is assumed to follow scaling relations derived 
from the spherical collapse model. Temperature T and halo mass 
M are assumed to be related by 



k B T 



M 



1/3 



6.03keV \l0 ,5 M Q /h) ^ ' ^\il(z)) \ 178 1 ' U> 

(e.g. lEke et alJl9~96llNavarro et alll995l) . The temperature k B T = 
6.03 k eV for a cluster with M = 10 15 Mq//i has been adopted from 
iMathiesen & Evrarcl fe oQlh . The density parameter at redshift z is 
denoted by £l(z), and A c is the mean overdensity of a virialised 
sphere, 

A c = 9n 2 {l+a[Q(z)-l] + Q(zf} (15) 

with (»,/?) = (0.7076, 0.4403) for a flat cosmology lStoehJl999h . 
Assuming that the total number N e of thermal electrons within the 
cluster virial radius is proportional to the virial mass yields 



N, 



1+/h . M 

~ .IB . 

2 m p 



(16) 



where / H is the hydrogen fraction of the baryonic mass (/ H ~ 0.76) 
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and m p is the proton mass. From X-ray data of an ensemble of 
45 clusters, Mo hr et alj fl999l) derived f B = 0.075/T 3/2 . Tradition- 
ally, the number density of dark matter hal oes is des cribed by the 
Press-Schechter formalism (Press & Schechter 1974). The comov- 
ing Press-Schechter mass function can be written as 



n PS (M,z) 



x exp 



V27r»+(z)M 2 

1 / M 



I n\ M 



(«+3)/6 



2D\{z) \M, 



(«+3)/3 



(17) 



where M, and p are the nonlinear mass today and the mean back- 
ground density at the present epoch, and D + (z) is the linear growth 
factor of density perturbations, normalised to unity today, D + (0) = 
i.n ~ -1 denotes the effective expon ent of the dark matter power 
spectrum at the cluster scale. ISheth & Tormenlfl9 99) recently pro- 
posed a significa ntly imp roved analytic derivation of the mass func- 
tion while Jenkins et al. 1 2001) measured the mass function of dark 
matter haloes in numerical simulations and found a fitting formula 
very close to Sheth & Tormen's, however, being of slightly lower 
amplitude at high masses. Thus, the fitting formula found by Jenk- 
ins et al. was used in our study. 

The total Compton-y parameter per unit solid angle is given 

by 



J/ n = d\ d 2 6y(6) : 



(18) 



where d A is the angular-diameter distance to the cluster. The mean 
background level of SZ fluctuations is given by 



<y bg > fl (z,M ) 



r d Z (i + Z ) 3 r 

J dz J M 
I dM I 



dMn PS (M,z0a(M,z) 



dVJ/ n (M,z) 



d 2 N(M,z) 
dMdV ' 



(19) 



where dV is the cosmic volume per unit redshift and unit solid 
angle, « PS (M, z) is the mass function of collapsed halos 1171 . and 
ifa(M, z) is the integrated Compton-y parameter per unit solid an- 
gle from 1181 expressed in terms of halo mass M and redshift z. 
Background fluctuations are due to Poisson fluctuations in the num- 
ber of clusters per unit mass and volume if cluster correlations are 
neglected. The variance of the background fluctuations reads 



( dM ( 

Jm J 



dV Wa{M,z)]' 



d 2 N(M, z) 
dMdV 



(20) 



Fig. shows a qualitative comparison of the influence of the 
background of unresolved SZ-clusters on the mean background 
level (ybg)s(z, M ) and the variance (y^)o(z, M ). For studying this 
influence quantitatively, the ratio of mean background levels and 
variances is defined via: 



(ybg>fl(z S im,^ 5 im) 
Cybg>e(Zraax,M m m)' 

(yt^eiz^m, M sim ) 

O^gMZmax, M mm )' 



(21) 



(22) 



where the numerator accounts for the resolved clusters in our sim- 
ulation with z < z s i m = 1.5 and M > M s]m = 5 x 10 13 Mq//j 
while the denominator accounts for all collapsed halos contribut- 



ing to the SZ-flux in our analytic estimate (z < Zn 



20 and 



M > M min = 1 x 10 13 Mq//;). Performing these integrals yields 
ratios of r mc;m = 40.6% and r v;lr = 93.3% and thus confirms the 
qualitative picture of Fig. [3] Therefore, we conclude that we can 




Figure 3. Comparison of the mean background level (yb g )e(z, 5 x 
10 13 M©//j) (solid) and the variances (y 2 } e (z, 5 X 10 13 M Q /h) (dashed) and 
(yl )e(z, \Q n M Q lh) (dash-dotted). The differential curves show qualita- 
tively the smaller impact of low-mass and high-redshift clusters on the vari- 
ance compared to the mean background of SZ fluctuations. 



safely neglect the effect of the background of unresolved SZ clus- 
ters on power spectrum statistics of our SZ all-sky map, especially 
when considering the mentioned uncertainties in / B and the ionisa- 
tion fraction of electrons in low-mass halos. 



5 RESULTS 

This section provides various characterisations of the SZ-cluster 
sample and properties of the resulting map. First, a visual impres- 
sion of the SZ-maps is given in Sect. 15,11 Distribution of the angular 
sizes and of the integrated thermal and kinetic Comptonisations are 
presented in Sect. l5.2l and in Sect. 15,31 respectively. The distribution 
of pixel amplitudes and a discussion of the sky-averaged thermal 
Comptonisation is given in Sect. 15 .41 The angular power-spectra in 
compariso n to those o btained in high-resolution simulations per- 
formed by White et al. 1 2002) is shown in Sect. 15.51 Finally, source 
counts in three relevant Planck-channels are given in Sect. 15.61 

In order to quantify the deviations resulting in using template 
SZ-maps instead of relying solely on analytical profiles and ide- 
alised scaling relations, the distributions following from the respec- 
tive approach are constrasted in Sect. l5.2l (angular sizes), Sect. 15.31 
(integrated Comptonisations) and Sect. 15.61 (source count at three 
selected P/ancA:-frequencies). 



5.1 Sky views 

In order to give a visual impression of the sky maps, all-sky views 
in Mollweide projection of the Compton-y parameter (Fig.0 as 
well as of the Compton-vv parameter (Fig. [5} are presented. Apart 
from those images, detailed maps of small regions of the SZ-sky 
are presented in Fig. [6] for the thermal and in Fig. [7] for the ki- 
netic SZ-effects, respectively. These detailed maps display interest- 
ing features: Clearly, cluster substructure is visible in the maps, e.g. 
at position (A,J3) ^ (134° 60, 45°25). 

Secondly, massive clusters that generate a strong thermal sig- 
nal, are rare, such that in drawing a peculiar velocity from a 
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Figure 4. All-sky map of the thermal Comptonisation parameter y in Mollweide projection. The shading is proportional to arsin h(10 6 Xy). 



Figure 5. All-sky map of the kinetic Comptonisation parameter w in Mollweide projection. The shading is proportional to arsinh(10 7 X w). 



Gaussian distribution large values are less likely to be obtained. 
Consequently, these clusters commonly show only a weak ki- 
netic signal, a nice example can be found at the position (A,/3) - 
(135°25,44°75). Closeby, the inverse example can be found at 
- (135° 40, 44°50), where a low-mass cluster shows only a 
weak thermal signal, but has sufficient optical depth and a high 



enough peculiar velocity to give rise to a strong kinetic signal. Fi- 
nally, at (A,{3) =! (135°90, 45°75), there is an example of a merging 
cluster, with a dipolar variation of the subcluster velocities. 

The occurence of high kinetic SZ-amplitudes is a subtle point: 
Cluster velocities follow a Gaussian distribution with mean con- 
sistent with zero, because the large scale structure is at rest in 
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Figure 6. Detail of the thermal Comptonisation map: A 2° X 2° wide cut- 
out centered on the ecliptic coordinates (A,/3) = (135°, 45°) is shown. The 
smoothing imposed was a Gaussian kernel with A0 = 2'.0 (FWHM). The 
shading indicates the value of the thermal Comptonisation y and is propor- 
tional to arsinh(10 6 xy). This map resulted from a projection on a Cartesian 
grid with mesh size ~ 14", i.e. no HEALPix pixelisation can be seen. 




13.6 13.8 14 14.2 14.4 14.6 14.8 15 15.2 15.4 15.6 

cluster mass log(M/(M //i)) 

Figure 8. Mean value and variance of the (Gaussian) distribution of peculiar 
velocities v r as a function of cluster mass M. The parameters of the distri- 
bution do not depend on mass, i.e. the mean is consistent with zero and the 
standard deviation has values - 320 km/s irrespective of mass. Standard 
deviations for the five bins corresponding to the largest cluster masses have 
been omitted due to poor statistics. The underlying data points represent 1% 
randomly selected entries of the Hubble-volume catalogue. 




ecliptic longitude .1 [degf 



Figure 7. Detail of the kinetic Comptonisation map: A 2° X 2° wide cut- 
out centered on the same position as Fig.|^| i.e. at the ecliptic coordinates 
(A,fS) = (135°, 45°) is shown. The smoothing imposed was a Gaussian ker- 
nel with A9 = 2'0 (FWHM). The kinetic Comptonisation w is indicated by 
the shading which is proportional to arsinh(10 X w). 



the comoving CMB-frame and with a standard deviation of o~ v = 
312.8 ± 0.2 km/s. This value has been measured for clusters in 
the Hubble-volume catalogue a nd is noticably smaller compared 
to <t v = 400 km/s proposed by Agh anim eta?] \ 1997ft . As Fig. [F] 
illustrates, the velocity-distribution does not depend on the clus- 
ter mass, because on the scales of typical cluster separation, linear 
structure formation is responsible for accelerating the clusters to 
their peculiar velocity. Massive clusters are rare and thus a high 
peculiar velocity is seldomly drawn from the underlying Gaussian 
distribution. Despite the seemingly large separation, it would be 
incorrect to draw the velocities independently from a Gaussian dis- 



tribution. Instead, the kinetic SZ-map ensures the consistency that 
the density and velocity fields have grown from the initial Gaus- 
sian random field by linear structure formation and have the correct 
relative phases. 

Similar to the clustering on large angular scales that the ther- 
mal SZ-map shows due to the formation of superclusters, the ki- 
netic SZ-map is expected to exhibit clustering on the same angular 
scales. This is because in the formation of superclusters, the veloc- 
ity vectors of infalling clusters point at the dynamical centre and 
are thus correlated despite the large separation. 

5.2 Distribution of angular sizes 

The distribution of cluster sizes is an important characteristic of the 
sky maps. For the derivation of core sizes, two different paths have 
been pursued in order to contrast the ideal case, in which cluster 
sizes follow from the well-known virial relations to the simulated 
and realistic case, in which the sizes are measured on the template 
maps themselves. First, t he cluster sizes are measured on th e data 
by fitting a King-profile ICavaliere & Fusco-Femianol 19781) to the 
thermal and kinetic template maps: 



y (r) = y • 



w (r) = Wo 



1 + 



(23) 



(24) 



yielding the core radii rj. f° r the thermal and rf for the kinetic 
map, respectively. The /3-flts have been centered on the pixel with 
the highest amplitude, and as free parameters only the central am- 
plitudes yo and wq were used apart from the core radii. The resulting 
radii have been averaged over all three projections of the cluster. 
Together with the comoving distance of the cluster as given by the 
Hubble-volume catalogue and the scaling factor required to match 
the size (compare Sect. l4.2> . the core radii have been converted into 
angular diameters. 
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Figure 9. The number of clusters n(O)d0 per square degree for given angular 
diameter 9 is shown without taking beam convolution into account, for ther- 
mal clusters (circles, dashed line) and kinetic clusters (crosses, dash-dotted 
line) as following from yS-model fits to template data. For comparison, the 
distribution of angular diameters obtained via the virial theorem (solid line) 
is also plotted. 



Secondly, an angular extent has been derived from the virial 
radius. Template data suggests the relation 



■ 0.12 r v : 



(25) 



rather than the value of r r ^ 0.07 r v ; r advocated by 
Llovd-Davies et alJ <2000l) and lKavetai] l200lh . In analogy, the 
angular diameter was then determined with the cluster distance 
given by the Hubble- volume catalogue. 

In Fig. 03 the size distributions for the thermal as well as for 
the kinetic clusters are given. Clearly, most clusters have angular 
diameters small compared to Planck's beam, and would appear as 
point sources. Here, it should be emphasised, that the HEALPix tes- 
selation with the chosen N Me -parameter does not resolve structures 
smaller than I'll, la the process of smoothing the clusters imposed 
prior to projection (compare Sect. I4.3> . clusters with diameters 
smaller than 2!0 have been replaced by 2-dimensional Gaussians 
with cr = 2'.0. Their normalisation corresponds to the integrated 
Comptonisations J/ and IV measured on the template maps. The 
smoothing is an absolute necessity because otherwise the HEALPix 
map would need ~ 10 4 times as many pixels for supporting the 
most distant and hence smallest clusters in the Hubble sample and 
hence ~ 10 4 times the storage space. A futher point to notice is the 
remarkably good agreement between diameters derived from the 
various prescriptions. 



5.3 Distribution of the integrated thermal and kinetic 
Comptonisation 

The signal strength of a cluster in an SZ observation is not given by 
the line-of-sight Comptonisation, but rather the Comptonisation in- 
tegrated over the solid angle subtended by the cluster. These quan- 
tities are refered to as the integrated thermal Comptonisation J/ and 
kinetic Comptonisation IV and are defined as: 



J/ 



dOy(0) and W 



I 



>>10" 




b o© 



integrated thermal Comptonisation log [if/ arcmin j 

Figure 10. The number of clusters «(J/)dJ/ per square degree for given 
integrated thermal Comptonisation J/ derived from template data (dashed 
line) in comparison to the analogous quantitiy based on virial estimates 
(solid line). 



For a simple model of the integrated Comptonisations as func- 
tions of cluster mass M, distance z and peculiar velocity u r it is as- 
sumed that the SZ-flux originates from inside a sphere of radius r vir , 
that the baryon fraction is equal to its universal value f B = Q. B /Q. M , 
that the ICM is completely ionised and has a uniform temperature 
predicted by the spherical collapse model laid down in eqn. 1141 . 
In this model, the actual distribution of electrons inside the virial 
sphere is of no importance. Then, the integrated Comptonisations 
are approximated by: 



3/\ 



1.98 v 
h 



'W-A 



0.29 



h 



Mi, 

Mvir 



d A \ ,, J Do 

Z) (l+z) o 



178, 



,(27) 



(28) 



dOw(0). 



(26) 



respectively. The reference values have been chosen to be Af» = 
10 15 M Q /h, d* = 100 Mpc/h and = 1000 km/s. d A is the an- 
gular diameter distance to the cluster. SI = Q(z) denotes the mass 
density at redshift z and A c = A c (z) the overdensity of a virialised 
sphere, an approximate description is given by eqn. H5\ . For typ- 
ical values for mass, distance and velocity, the thermal and kinetic 
SZ-effects differ by approximately one order of magnitude. The 
baryon fraction is set to the universal value f B = Q. B /Q M = 0.133 
for the remainder of this paper. 

Distributions of the integrated thermal and kinetic Comptoni- 
sations are shown in Figs. llOl and ll II respectively. The distributions 
have been derived from actual scaled template data in comparison 
to the values obtained from 1271 and 12 SI . 

Fig. llOl shows the number of clusters per deg 2 with integrated 
thermal Comptonisation J/. It can be seen that the approach rely- 
ing on the virial theorem underestimates the number of cluster by a 
factor 2-3 for large integrated Comptonisations. Alternatively, one 
could state that the distributions are separated at high Comptonisa- 
tions by slightly less than 0.5 dex. The reason for the significantly 
larger integrated Comptonisations determined from template data 
is due to the fact that the template clusters were matched to the 
catalogue entries given by the Hubble-volume simulation by their 
virial masses. In irregular clusters, there is a significant fraction 
of the gas located outside the virial sphere and thus the integrated 
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integrated kinetic Coraptonisation log I TV/ arcmin J 

Figure 11. The number of clusters n( r W)(¥W per square degree for given 
integrated kinetic Comptonisation *W derived from template data (dashed 
line) in comparison to the analogous quantitiy based on virial estimates 
(solid line). Here, the scaling of the vertical axis is linear, in contrast to 
Fig. II II such that the underlying Gaussian distribution of peculiar velocities 
becomes apparent. 
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Figure 12. Distribution of pixel amplitudes of the thermal Comptonisation 
parameter y (circles, solid) and of the absolute value of the kinetic Compton- 
isation | vv ■ | (squares, dashed). Zero values have been deliberately excluded. 
The figure illustrates the occurence of clusters comparable to the HEALPix 
pixel scale: The thin set of lines shows an additional peak at small Compton- 
isations, that vanish after convolution with a beam of A8 = 5'.0 (FWHM), 
i.e. comparable to Planck, as shown by the thick lines. 



Comptonisation is systematically underestimated when applying a 
spherical overdensity code to simulation data, as previously exam- 
ined in Fig. [2] In scaling the template clusters up to the masses 
required by the Hubble-volume catalogue, this difference is am- 
plified because in the sparsely sampled region of the M-z-diagram 
(compare Figs. Q and |2j clusters have on average to be scaled to 
higher masses, which explains the offset in the distributions. A sec- 
ond effect is the evolution of ICM temperature. Compared to the 
temperature model eqn. J 141 based on spherical collapse theory, 
the plasma temperatures are smaller by approximately 25%, i.e. 
the mean SPH-temperature of the particles inside the virial sphere 
is smaller than expected from spherical collapse theory and reflects 
the departure from isothermality: The template clusters do show a 
temperature profile that declines towards the outskirts of the clus- 
ters, which decreases the integrated Comptonisation relative to the 
values derived by means of the virial theorem. 

Furthermore, the dependence of electron temperature on clus- 
ter mass is noticably weaker than the M 2/3 -scaling: The cluster 
number weighted average for the exponent a in the scaling T oc M a 
relating temperature to mass was found to be (or) = 0.624, and at 
the redshifts around unity, where most of the clusters reside, val- 
ues as small as a = 0.605 were derived. Using this scaling, the 
Compton-y parameter and hence the integrated thermal Compton- 
isation if shows a significantly shallower distribution compared to 
the distribution relying on simple scaling arguments. 

The same argument applies to the kinetic Comptonisation "W, 
as depicted in Fig. ^3 Here a shifting of the values to smaller ki- 
netic Comptonisations is observed when comparing estimates fol- 
lowing from the virial theorem to actual simulation data. The shift 
of the peak of the distribution amounts to about one third dex, 
as explained above for the thermal Comptonisations. Keeping the 
-scaling of the thermal SZ-effect in mind, the shift in the TV- 
distribution is then consistent with the shift of the J/-distribution. 



5.4 Distribution of Comptonisation per pixel 

Fig. ll2l shows the distribution of the pixel amplitudes of the thermal 
SZ-map as well as of their absolute values in the kinetic SZ-map. 
Clearly, the kinetic and thermal SZ-effects are separated by approx- 
imately one order of magnitude. 

The distribution of pixel amplitudes is very broad, encompass- 
ing the largest line-of-sight Comptonisations of y - 1.5 x 10~ 4 and 
I vv| =! 1.6 x 10~ 5 down to very low signals below logfj) m -20. The 
distribution is bimodal, which is a pixelisation artefact and which 
is caused by the replacement of faint and small clusters with a very 
narrow Gaussian, the extent of which is slightly above the pixel 
scale, once the cluster is smaller than = 2!0 in diameter. These 
clusters are more concentrated than the King-profiles of resolved 
clusters. There is a caveat when applying an expansion into spheri- 
cal harmonics Y c m {0) to the SZ-maps: The smallest clusters are only 
a few pixels in diameter. Working with the HEALPix tesselation, 
reliable expansion coefficients can only be obtained up to multi- 
pole moments of order t 2 x A'side, i.e. up to I = 4096 in our case, 
which corresponds to angular scales of 2'64. The pixel scale is 1 '7 1 
for this choice of the -parameter. Consequently, the small clus- 
ters will not be contained in an expansion into spherical harmonics, 
as shown by the set of thick lines in Fig. 1121 Here, the map has 
been decomposed in ai m -coefficients (compare eqn. <31» . multi- 
plied with the fl«)-coefficients of a Gaussian beam of A8 = 5'.0 
(FWHM) and synthesised again. Then, the resulting smoothed map 
does not contain small clusters, because the decomposition into 
spherical harmonics has not been able to resolve structures that ex- 
tend over only a few pixels. 

The mean value of the thermal Comptonisation y has been 
determined to be (y) = 3.01 x 10~ 7 and the pixel-to-pixel vari- 
ance is cr y = -\/{y 2 ) - (y) 2 = 1.85 x 10~ 6 . In analogy, the value 
(w) = 6.28 x 10~ 9 has been derived for the kinetic map, with vari- 
ance <x„, = -\/{w 2 ) - (w) 1 = 3.78 x 10~ 7 , i.e. the mean kinetic 
Comptonisation is consistent with zero, due to the peculiar veloc- 
ities following a Gaussian distribution with zero mean. The mean 



All-sky thermal and kinetic SZ-maps for Planck 1 1 



value of the moduli of the pixel amplitudes of the kinetic map is 
<|vv|> = 7.65x 1CT 8 . 

The value for the mean Comptonisation (y) measured on the 
map should account for roughly 40% of the mean thermal Comp- 
tonisation as derived in Sect. 14.41 due to the lower mass thresh- 
old inherent to the simulation. Keeping in mind the absence of 
any diffuse component of the thermal Comptonisation, the value 
deriv ed here is compatible with the value of ^ 10~ 6 given by 
iRefregier et alJ fcOOd) and [ Scaramella"et"aT] ll993l) . but falls short 
of the value derived by White et al. (2002) by a factor of less than 
two. lMvers et alJ 120041) performed a cross-correlation of WMAP- 
data with clusters from the APM survey and found the mean 
Comptonis ation to be significantly larger and to be in accordance 
with lPersi et alJ 11995). but in contradiction with expectations from 
CDM models. 




5.5 Angular power spectra of the thermal and kinetic 
SZ-effects 

In this section, the angular power spectra are given for the all-sky 
maps. They follow from a decomposition of the spherical data set 
into spherical harmonics Y'"(Q): 
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(29) 
(30) 
(31) 



respectively. The spherical harmonical transform w'. has been de- 
termined from the absolute values of the kinetic map amplitudes. 
The reason for doing so is the vanishing expectation value of the 
peculiar velocities in the comoving frame such that for a given clus- 
ter in the thermal SZ-map, both signs of the kinetic SZ-effect are 
equally likely to occur and the cross-power averages out to zero. 
The angular power spectra and the cross power spectrum are de- 
fined via: 
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(32) 
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(34) 



with the asterisk denoting complex conjugation. The resulting 
power spectra are g iven in Fig. 1131 in comparison to the power- 
spectra derived by White et al. (20Q2) in simulations covering 
smaller angular scales. The curves match well, and the remaining 
discrepancies may be explained by the fact that in the maps pre- 
sented here, power is missing on small scales due to the low-mass 
cutoff, whereas the simulation by White is missing power on large 
scales due to the smallness of their simulation box. The bending- 
over of the spectra derived from our SZ-maps is also due to the fact 
that the expansion in spherical harmonics cannot be computed for 
angular scales approaching the pixel scale and thus does not include 
very small clusters of sizes comparable to the pixel size, as already 
discussed in Sect. 15.41 

If clusters were randomly positioned point-sources on the sky, 
the number of clusters per solid angle element would be a Poisson- 
process and the resulting power spectrum should be flat, i.e. C(£) oc 



multipole order £ 

Figure 13. Angular power spectra of the thermal and kinetic SZ-effect: 
( ', , (7 1 (circles, solid line), C ww (£) (crosses, solid line) and the cross power 
spectrum C yw i(f) (diamonds, solid line) are shown in comparison to the 
power spectra of the thermal SZ-effect (circles, dashed line) and the kinetic 
SZ-effect (crosses, dashed line) obtained by White et al. 1 2002) at smaller 
scales, i.e. at higher multipole order t. 



N (N ist the number of sources), as shown bv lSeliak et alJ 1200 lh . 
Contrarily, the brightness distribution of clusters assigns additional 
weight to the large angular scales and giving rise to a significant 
deviation in the slope of the power spectra C{£) oc £~ y as a function 
of £: The measured slope is 2 + y = 1.53 ± 0.07 for the thermal and 
2 + y = 1.45 ± 0.07 for the kinetic SZ-effect, which reflects the de- 
viation from pure Poissonianity. In the fitting, the values for £ have 
been restricted to 1 < £ < 100 and the errors derived correspond to 
the 95% confidence intervals. 

Furthermore, Fig. ll3l shows the cross-correlation between the 
thermal SZ-map and the absolute value of the kinetic SZ-map. As 
expected, the amplitude of the cross-power spectrum is at an in- 
termediate level compared to autocorrelations of the thermal and 
kinetic SZ-maps. 



5.6 Source counts at Planck frequencies 

As the last point in this analysis, we address the SZ source counts, 
i.e. the number N of SZ-clusters giving rise to flux changes ex- 
ceeding a certain flux threshold S m i„. The SZ flux modulation as a 
function of frequency is given by: 



S(x) 



SaJdSl [yg(x)-pTh(x)] 

So [J/gW - 'Wh(x)] = Sy(x) - sve*), 



(35) 
(36) 



where Sq = 22.9 Jy/arcmin~ is the flux density of the CMB and 
if and "W denote the integrated thermal and kinetic Comptonisa- 
tions. The functions g(x) and h(x) are the flux modulations caused 
by the thermal and kinetic SZ-effects for non-relativistic electron 
velocities: 



g(x) 
h(x) 



x 4 exp(x) 
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Figure 14. Source counts N((S) > S m \ n ) for thermal (circles) and kinetic 
clusters (crosses) for Planck's vq = 143 GHz channel and for fluxed mea- 
sured on the scaled template clusters (dashed line) in comparison to virial 
fluxes (solid line). 



Here, x again denotes the dimensionless frequency x = 
hv/(k B T C MB)- The averaged flux (S) vo at the fiducial frequency v 
is obtained by weighted summation with the frequency response 
window function R vo (v) and can readily be converted to antenna 
temperature T A by means of eqn. I39i : 



<S>v„ 



J dvS(v)R V0 (v) 
/ dvR V0 (v) 



2-\k B T A 



(39) 



The main characteristics of Planck's receivers and the conver- 
sion factors from 1 arcmin 2 of thermal or kinetic Comptonisation 
to fluxes in Jansky and changes in antenna temperature measured 
in nK is given by Table|5| For the derivation of the values a top-hat 
shaped frequency response function R vo (v) has been assumed: 



Rv (v) ■ 



1, V 6 [v 
0, V <t [V 



Ay, v + Ay] 
Ay, y + Ay] 



(40) 



Figures fT4ll 15l and l 1 6l show the source counts stated in number 
of clusters per deg 2 as a function of averaged flux (S > vo for Planck's 
v = 143 GHz-, v = 217 GHz- and y = 353 GHz-channels, 
respectively. 

The source counts N(S ) are well approximated by power laws 
of the form: 



N((S)>S mln ) = N S". 



(41) 



Values for the normalisations No and the slopes a have been ob- 
tained by fits to the source counts for the three relevant Planck- 
frequencies and are stated in Tabled In the fits, the four rightmost 
bins have been excluded because of poor statistics. The parameters 
of the power law has been derived for the fluxes following from the 
idealised case based on the virial theorem and compared to fluxes 
determined from template cluster data. 

The slopes derived from fits to the cluster number counts are 
slightly steeper for the virial estimates compared to template data 
(a a -5/3 versus a a -1.4), which again reflects the weaker de- 
pendence on cluster mass observed in template data. Comparing 
data sets for different frequencies, the slopes a are of course almost 
identical, because only amplitudes are changed by the choice of a 
different frequency band. The number of clusters No stays roughly 
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Figure 15. Source counts N((S) > S mm ) for thermal (circles) and kinetic 
clusters (crosses) for Planck's vo = 217 GHz channel, the dashed and solid 
lines contrast the fluxes measured on the template data and those follow- 
ing from virial scaling relations, respectively. For the given frequency re- 
sponse function K,, (v), the thermal SZ-effect does not vanish entirely at 
v = 217 GHz. 
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Figure 16. Source counts N({S) > 5 m ; n ) for thermal (circles) and kinetic 
clusters (crosses) for Planck's vo = 353 GHz channel, again for fluxes de- 
rived from template data (dashed line) in comparison to fluxes following 
from virial scaling relations (solid line). 



constant in the case of the kinetic SZ-effect, but reflects the dis- 
tinct frequency modulation in the case of the thermal SZ-effect. 
Here, it should be emphasised, that the thermal SZ-effect does not 
vanish entirely at v = 217 GHz due to Planck's symmetric fre- 
quency response window. The difference in numbers between the 
estimates based on virial quantities to those measured on template 
data amounts to roughly half an order of magnitude in the kinetic 
SZ-effect, but rises almost an order of magnitude at small fluxes 
for the thermal SZ-effect. There is however good agreement in the 
number counts of thermal SZ clusters at high fluxes. 

The difference in slope of the thermal versus kinetic SZ- 
cluster number counts is caused by the M 5/3 -scaling of the thermal 
SZ-effect relative to the proportionality to M of the kinetic effect. 
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Planck channel 


1 


2 


3 


4 


5 


6 


7 


8 


9 


centre frequency vo 


30 GHz 


44 GHz 


70 GHz 


100 GHz 


143 GHz 


217 GHz 


353 GHz 


545 GHz 


857 GHz 


frequency window Ay 


3.0 GHz 


4.4 GHz 


7.0 GHz 


16.7 GHz 


23.8 GHz 


36.2 GHz 


58.8 GHz 


90.7 GHz 


142.8 GHz 


thermal SZ flux (Sy) 


-12.2 Jy 


-24.8 Jy 


-53.6 Jy 


-82.1 Jy 


-88.8 Jy 


-0.7 Jy 


146.0 Jy 


76.8 Jy 


5.4 Jy 


kinetic SZ flux (S^y) 


6.2 Jy 


13.1 Jy 


30.6 Jy 


55.0 Jy 


86.9 Jy 


110.0 Jy 


69.1 Jy 


15.0 Jy 


0.5 Jy 


antenna temperature ATy 


-440 nK 


-417 nK 


-356 nK 


-267 nK 


-141 nK 


-0.5 nK 


38 nK 


8.4 nK 


0.2 nK 


antenna temperature AT>w 


226 nK 


220 nK 


204 nK 


179 nK 


138 nK 


76 nK 


18 nK 


1.6 nK 


0.02 nK 



Table 2. Characteristics of Planck's LFI- and HFI- receivers: centre frequency v, frequency window Av (as defined in eqn. 1401 ). fluxes (Sy) and (Siy) (see 
eqn. <36l ) generated by the respective Comptonisation of J/ = *W = 1 arcmin 2 and the corresponding changes in antenna temperature ATy and AT-w- Due to 
Planck's symmetric frequency response window, the thermal SZ-effect does not vanish entirely at v = 217 GHz. 



Planck channel v = 143 GHz v = 217 GHz v = 353 GHz 



thermal SZ-effect, virial estimate 


lOE 


\N = 


-1.78 


± 


0.02 


lOE 


;N = 


-5.31 ± 


0.02 


lOE 


;M> = 


-1.42 


± 


0.22 






a = 


-1.67 


+ 


0.03 




a = 


-1.67 ± 


0.07 




a = 


-1.66 


+ 


0.02 


kinetic SZ-effect, virial estimate 


lOE 


\N = 


-4.49 


± 


0.01 


lOE 


\N = 


-4.42 ± 


0.01 


lOE 


\N = 


-4.57 


± 


0.01 






a = 


-0.76 


± 


0.06 




a - 


-0.76 ± 


0.05 




a = 


-0.76 


± 


0.06 


thermal SZ-effect, simulation 


lOE 


;N = 


-2.36 


+ 


0.02 


lOE 


;N = 


-5.31 ± 


0.02 


lOE 


;N = 


-2.06 


± 


0.02 






a = 


-1.40 


+ 


0.03 




a = 


-1.40 ± 


0.06 




a = 


-1.40 


+ 


0.03 


kinetic SZ-effect, simulation 


lof 


lN = 


-3.95 


+ 


0.01 


lOE 


lN = 


-3.88 ± 


0.01 


lOE 


lN = 


-4.02 


+ 


0.01 






a = 


-0.72 


+ 


0.04 




a = 


-0.72 ± 


0.05 




a = 


-0.72 


+ 


0.05 



Table 3. Values obtained from fits of a power law of the type N(S ) = No S a to the cumulative source counts as a function of flux exceeding the threshold S for 
both SZ-eft'ects. In the table, values obtained from virial estimates are contrasted to values following from measurements on template data. The errors quoted 
denote the 95% confidence intervals. 

within the boundaries in mass (M > 5 x 10 13 Mq/K) and red- 
shift (z < 1.48) make up = 40% of the mean Comptonisation, 
but account for ~ 98% of the variance. The value for the mean 
Comptonisation corresponds well to that obtained by other authors 
(Sect.l5~4l 

• The power spectra (Sect .l5.5t are compatib le in amplitude and 
slope to the ones found bv lWhite et alj J2002I) . On large angular 
scales, i.e. at small multipole I, deviations from the Poissonianity 
in the slope of the power spectrum have been found. 

• The velocities of the kinetic SZ-map correspond to the actual 
cosmological density environment, i.e. correlated infall velocities 
are observed due to the formation of superclusters, which highlights 
a significant improvement in comparison to methods that draw a 
cluster peculiar velocity from a (Gaussian) distribution and enables 
searches for the kinetic SZ-effect by considering spatial correla- 
tions with the thermal SZ-effect. The cross correlation of the ther- 
mal with the kinetic SZ-map yields a spectrum similar in shape at 
intermediate amplitudes (see Sect. 15. 5> . 

Despite the high level of authenticity that the all-sky SZ-maps 
exhibit, there are quite a few issues not being taken account of: The 
baryon distribution and temperature inside the ICM is governed by 
processes beyond adiabatic gas physics, for example in the form 
of supernova feedback and radiative cooling. Especially the lat- 
ter process gives rise to cool cores which may enhance the ther- 
mal SZ-signal. The ionisation inside the clusters was assumed to 
be complete. Furthermore, the maps contain only collapsed objects 
and hence filamentary structures or diffuse gas are not included. 
Concerning the thermal history of the ICM, reionisation had to be 
neglected. The kinetic map has been constructed without taking ac- 
count of velocity fluctuations inside the cluster. This does not pose 
a problem for Planck, but needs to be remedied in high-resolution 
SZ-surveys to be undertaken with the Atacama Cosmology Tele- 



Due to the difference in slope, the effects are separated by two or- 
ders of magnitude for the largest fluxes, while this difference in- 
creases to eight orders of magnitude for the smallest fluxes, which 
hints at the difficulties to be faced in detecting kinetic clusters com- 
pared to even faint thermal detections. The slopes derive d here are 
in good agreement with the those obtained by Kav et al. 1 2001) 



6 SUMMARY 

All-sky maps for the thermal and kinetic Sunyaev-Zel'dovich ef- 
fects are presented and their characteristics are described in detail. 
The maps because of their angular resolution and the data stor- 
age format chosen (HEALPix) especially suited for simulations for 
Planck. 

• The all-sky maps of the thermal and kinetic SZ-effects pre- 
sented here incorporate the correct 2-point correlation function, the 
evolving mass function and the correct size distribution of clusters, 
to within the accuracy of the underlying Hubble-volume simulation 
and the small-scale adiabatic gas simulations. 

• The maps presented here exhibit significant cluster substruc- 
ture (compare Sect. I5.lt . In spite of this, fits to the Comptonisation 
maps yield angular core radii, the distribution of which are close to 
the expectation based on the virial theorem (Sect. l5~2l . 

• The difference in the distribution of the integrated Compton- 
isations J/ and "W (Sect. I5~3l and source counts N((S) > S m in) 
(Sect. l5.6l between values derived from scaling relations compared 
to those following from template data have been found to be sub- 
stantial, which hints at possible misestimations of the number of 
clusters detectable for Planck. 

• An analytic investigation in Sect. l4.4| quantified the contribu- 
tion of the cluster sample to the sky averaged mean thermal Comp- 
tonisation (y) and its variance <r v . It was found that the clusters 
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scope and the South Pole Telescope. Yet another imperfection is 
the lack of non-therm al particle populations that cause the rela- 
tivistic SZ-effect I EnBlin & Kais eJ2000h . whose detectability with 
Planck is still a matter of debate. 
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